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Abstract. We investigate a non-locality of Moss-Okninski transforma- 
tion (MOT) used to separate positive and negative energy states in the 3+1 Dirac 
equation for relativistic electrons in the presence of a magnetic field. Properties of 
functional kernels generated by the MOT are analyzed and kernel non-localities 
are characterized by calculating their second moments parallel and perpendicu- 
lar to the magnetic field. Transformed functions are described and investigated 
by computing their variances. It is shown that the non-locality of the energy- 
separating transformation in the direction parallel to the magnetic field is char- 
acterized by the Compton wavelength Ac = h/mc. In the plane transverse to 
magnetic field the non-locality depends both on magnetic radius L = (h/eB)^^^ 
and Ac. The non-locality of MO transformation for the 2+1 Dirac equation is also 
considered. 
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1. Introduction 

The Dirac equation, in spite of its numerous applications and fundamental significance 
in relativity and quantum mechanics, is far from being completely understood. This 
is especially true for the Dirac electrons in the presence of external fields. For 
example, it has been known for many years that the Dirac equation can be solved 
exactly in case of an external uniform magnetic field [1], but it is only recently that 
an explicit electron dynamics of Dirac electrons in the presence of a magnetic field 
was worked out [2]. This problem is related to a somewhat mysterious phenomenon 
known in the literature under the name of Zitterbewegung. The phenomenon was 
conceived in 1930 by Schrodinger who observed that, according to the Dirac equation, 
operators corresponding to the electron velocity do not commute with the Dirac 
Hamiltonian, so the velocity depends on time also in absence of external fields [3]. 
It is understood at present that this remarkable phenomenon is due to an interference 
of states corresponding to positive and negative electron energies [31 [6]. In fact, it 
is known by now that, whenever one deals with a spectrum of positive and negative 
energy branches, the Zitterbewegung will appear [7|. 

There exist ways to avoid the duality of positive and negative energies. The best 
known of them is the Foldy-Wouthuysen transformation (FWT) which, in case of no 
external fields, allows one to break exactly the 4x4 Dirac equation into two 2x2 
equations for positive or negative energies [S]. In their pioneering paper, FW remarked 
that a kernel transforming functions from the original to the FW representation is 
characterized by a non-locality in coordinate space on the order of the Compton 
wavelength Ac = h/mc. Rose in his book W put this statement on a quantitative 
basis by showing that the second moment of the functional kernel is equal to (3/4)A^. 
This result was not followed by other investigations and it is by now not well known. 
The problem was recently taken up by the present authors who demonstrated that also 
other transformed quantities are characterized by non-localities on the order of Ac [10] . 

Case 11 showed that also the Dirac equation including the presence of an external 
uniform magnetic field B can be exactly separated into two equations for positive and 
negative electron energies. This form allows one to find easily the energy eigenvalues 
in any desired semi-relativistic approximation (see our Eq. ([7]) and Ref. [51]). It is 
possible to separate the positive and negative electron energies for i? 7^ because the 
presence of a uniform magnetic field does not induce electron transitions between the 
two parts of the spectrum. The case of magnetic field is of interest for several reasons. 
First, it introduces an adjustable external parameter B affecting all energies and other 
electron characteristics. Second, it quantizes the orbital motion in two degrees of 
freedom. Third, it quantizes the spin energies. However, to our knowledge properties 
of the separating transformations for this important case were not investigated until 
now. 

In the form proposed by Case the transformation operator is C/ = e^^^'^/'^ with p2 
being the standard Dirac matrix, and the phase given hy (f) — tan~^(o'7r)/mc, where & 
are the Pauli matrices and tt is the generalized momentum operator including the 
vector potential of the magnetic field. In this form it is difficult to analyze the 
properties of U . A method to get a more explicit form of U was described by de 
Vries [12]. However, it was shown that the "separating" transformations, both in 
absence of fields and for B 7^ 0, are not unique. In other words, it is possible to 
separate the Dirac equation into two equations for positive or negative energies, 
both in absence of fields and in the presence of a magnetic field, using different 



Non-locality of transformations in a magnetic field 



3 



unitary transformations [T^ [T51 UM [T5] . In our description involving Dirac 
electrons in a magnetic field we use a two-step transformation proposed by Moss and 
Okninski (MO) [H]. An advantage of this choice is that the MOT is given explicitly, 
it is convenient for kernel calculations and we can compare our final results in the 
limit i? — > with our previous description for _B = 0, see |10| . In our treatment we 
do not concentrate so much on the separation of energies but rather on the properties 
of transformation operators and the transformed functions. Naturally, we emphasize 
the effects of an external magnetic field. 

Our subject is of relevance for at least three reasons. First, with the comeback of 
interest in narrow-gap semiconductors and, in particular, the appearance of zero-gap 
graphene [16], there has been a real surge of works concerned with relativistic-type 
wave equations. Second, it is now possible to simulate the Dirac equation with the 
trapped ions or cold atoms interacting with laser radiation, where one can tailor 
"user friendly" values of the basic parameters mc^ and c [H [T71 [TSJ [H]. In fact, a 
proof-of-principle experiment simulating the 1-1-1 Dirac equation and the resulting 
Zitterbewegung was carried out by Gerritsma et al. [20]. Third, one can simulate 
equally well the Dirac equation in the presence of a magnetic field [21 [21] . In such 
simulations it is possible to tailor a seeming intensity of the field making its effects on 
relativistic electrons relevant in terrestrial conditions. Putting it in more quantitative 
terms, one can achieve experimentally the so-called Schwinger field, defined by the 
equality heBo/m = mc^, not only in the vicinity of neutron stars but in standard 
laboratories, see our discussion below. 

Our paper is organized in the following way. In Sections II and III we introduce the 
Moss-Okninski transformation, describe its functional kernel and calculate its zeroth 
and second moments. In Section IV we transform Gaussian functions and calculate 
their normalized variances. In Section V we analyze kernels in two dimensions. In 
Section VI we discuss our results. The paper is concluded by a summary. In two 
Appendices we give some mathematical details of our derivations. 



2. Transformation Kernel and its non-locality 



In this section we define the MO transformation and its functional kernel. Then we 
analyze kernel's properties and, in particular, its non-locality. The Dirac Hamiltonian 
for a relativistic electron in a magnetic field is 

"H ~ COixTTx + CCtyTTy + CCtzT^z + ^mC^ , (1) 

where tt = p + 2^ is the generalized momentum, |e| is the electron charge, A is the 
vector potential of a magnetic field B, and di and /3 are the Dirac matrices in the 
standard notation. Following Moss and Okninski [15) we first introduce a unitary 
operator 

(2) 



in which 
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This operator transforms the initial Hamihonian into a form with vanishing diagonal 
terms 

= caxTTx + cayTTy + ca^TTz + Smc . (4) 
The second step transforms "H to a strictly diagonal form with the use of an operator U 

?7 = (/3 + e) , (5) 

where 

e = „^ , . (6) 

(H'2)l/2 

After the complete transformation the Hamiltonian has a diagonal form (cf. Refs.|lll 

UnU^ = Pim^c^ + c^TT^ - ec^aSB)!/^, (7) 
where S = ( ^ ) is the 4x4 spin operator. Since tt is an operator, the square 



X cr ^ 

root is understood as an expansion in powers of tt. In the following we concentrate 
on the properties of the transformation and transformed functions. 

Let l^") be a function in the old space and \^') — V]"^) in the transformed space. 
We introduce standard projecting operators for positive and negative energy states 

F± = i(i±e). (8) 

If the projection operators are applied to the transformed functions l^*'), they split 
them into states corresponding to positive or negative energies: |^''^) — P^\^'). 
Next l^*'^) are transformed using the operator t/, i.e. |^'^) — UY^'^). We introduce 
operators IJ^ defined as = JJP^ with the property 1^'=*=) = IJP^\^'). We have 

c>± = f /3 + e) f i ± e 



2^2 

Let us consider transformed functions \f'^(ri) = (ri|?7*|5''). Inserting the identity 
operator 1 = / \r2){r2\<P'r2 one has 

k^{rur2W{r2)d?r2, (10) 

where we defined the functional transformation kernels 

k^{ri,r2) = {ri\U^\r2). (11) 

In the following we investigate the properties of K^{ri,r2)- To calcu- 
late K'^{ri,r2) we insert twice the unity operator 1 = |m)(m|. The states |m) 
are eigenstates of the Hamiltonian Ti, given in Eq. ([4]) , and the summation is carried 
out over all quantum numbers characterizing |m). Thus 

i^±(ri,r2) - Yl (n|ni)(m|;7±|m')(m'|r2). (12) 
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The operator consists of operators: i, (5 and Q, see Eq. The matrix elements 
of /3 and 1 are diagonal: /3m, m' — PSm.m' and l(5m,m', respectively. The non-diagonal 
contribution to f{^{ri,r2) arises from 9 

e(ri,r2) - (ri|m)(m|e|m')(m'|r2). (13) 

m,m' 

The eigenstates |m) of the Hamiltonian H are related to the eigenstates |n) of the 
initial Dirac Hamiltonian H by a unitary transformation 

|m) = V\n), (14) 

where V is defined in Eq. ([2]). We want to calculate the elements of 0(T'i,r2) matrix 
and for this purpose we need to know explicit forms of the eigenstates |m). 

We take the vector potential of a magnetic field in the Landau gauge A = 
{—By, 0, 0). Then the eigenstates |n) depend on five quantum numbers: n, kx, kz, e, s, 
where n is the harmonic oscillator number, kx and are wave vector components, e = 
±1 labels positive and negative energy branches, and s — ±1 is the spin 
index corresponding to the spins ±1/2, respectively. In the Johnson-Lippman 
representation [22] the eigenfunctions ^ni^) = (''In) of the Hamiltonian H are 

f si{eEn^k, + mc^) (t>n-i{r) \ 
S2{eEn^k, + mc^) (/>ri(r) 

(siCPz - S2huJn) (/>«-i(t-) 
\ -{sihuJn+ S2CPz) (/>n(r) J 

where si = (1 + s)/2 and S2 = (1 — s)/2 select the states s = ±1, respectively. The 
frequency is aj„ = ujy/n with uj = \/2c/L, the energy is 



nepz 



(15) 



En,k^ = V(mc2)2 + (nw„)2 + (UkzC)^, (16) 
and it does not depend explicitly on s. The norm in Eq. ((T5|) is Nnek^ — i'^^n k + 
2emc^£'ri,fe^)~^/^, and the (/>„ functions are 

Mr)^ 2^7ZC ^-^^^^ ' 



where H„(^) are the Hermite polynomials, C„ = •\/2"n!^/7r, the magnetic radius 
is L = \fhJeB and ^ = yjL — kxL. 

Since V defined above is an unitary matrix, there is -Em — E^ — eEn^k^- The 
states |m) are characterized by the same five quantum numbers describing |n), so there 
is also (5m, m' = ^n,n'- For two eigcnstatcs |m) and |m') of the Hamiltonian T-L we have: 
{m\n'\m') = Em(5m,m', {m\H'^\m') = EfJ^^n,', and 

E 

(m|e|m') = (5m, m' = qeSnAi', (18) 

\/ E 



where q — 1 or q = —1. Selecting q = 1 we obtain from Eqs. (fT5)) . and (HI 

(nieira) =^e(ri|m)(m|r2) =^e(ri|\/n)(nf|r2). (19) 

m n 

Selecting = — 1 in Eq. would not have changed the results because we perform 
a summation over e — ±1. Since 1^ is a number matrix, there is (r|yn) = V{r\n), and 
one can rewrite (t'i|6|t2) in the form 

{ri\e\r2) =Veoiri,r2)V\ (20) 
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where we introduced an operator Qq by the relation 

(ri|eo|r2)=^6(ri|n)(nh). (21) 

n 

The advantage of 6o('"i, '"2): as compared to 6(ri,r2), is that 6o('"i,'''2) can be 
directly expressed in terms of the eigenstates |n) given by Eq. (ITS]) . Writing in Eq. (PT|) 
all quantum numbers we have 

/■oo /"OO 

(ri|eo|r2) = / dkj dfc. ^ e$n(n)$l;(r-2), (22) 

where $n(^) are given in Eq. (|15p . The integrals over and as well as the 
summations over e and s may be done analytically, the summation over n is to be 
performed numerically. However, in the limit of strong magnetic fields analytical 
results may be obtained. The summations and integrations indicated in Eq. ()22|) are 
carried out in the following order. First, we sum over the spin variable s noting 
that S1S2 = and J2s^i ~ 12s^2 — 1- I^i the summation over e we use the 
property: J2e=±i e = 0. As a result, one obtains 



poo I* 

eo(ri,r2) = ^ / dkJ 



dk. 



,(23) 



eo(ri,r2) = £ 



/ 'yi(f>n'{ri)<Pl,{r2) 73'/>n'(ri)<?!);,(r2) -720n' (r-i)*/),* (r2) \ 

7i(/)„(ri)(/),*(r2) ~720n(r-i)C'('^2) -730n(ri)0;(r2) 

73(f>n'{ri)<j)*^,{r2) -720n'(n)0;(r-2) -Ji(l)n'{ri)(f>*^,{r2) 

V -I2<l>n{ri)<l>n'{r2) -73 (n )(/),* (r2) -7i0n(ri)C(r2) / 

where 71 = mc^ /{2ttE), 72 = ?icij„/(27riJ), 73 = cpz/i^nE) and = n — 1. The 
matrix elements of 0o(t'i, r2) are products of integrals over k^ and k^. There appear 
three different integrals over k^: Di, D2 and D3 and four integrals over k^ denoted 
as Am,n- We write explicitly 

^ i:>lA„_i,„_l D^An-l^n-l -y/nD2An-l.n \ 

n,n — 

£'3A„-1,„-1 ~^D2An-\,n ~DxAn-\^n-\ 

\ -y/nD2An^n-\ -D-iAn^n -DiA„_„ / 

The integrals over k^ are 

Am.n= / e-M"^ 7^ H„,(a)H„(g2)e'^/'^^'+^^^dfc., (25) 
where = yj / L — k^L with j = 1,2. They are calculated analytically 

A„-,,„-, = ;^L^^.(p^)e-'^V2+.x^ (26) 

A„_i,„ =-^L^,(p-^)e-'^^/^+^^ (27) 

A„,„_i =^L^,(p-^)e-P^/2+»^ (28) 

where j = 0, 1, p = [(a;i — 2:2)^ + (yi — ?/2)^]/(2i^), L^(p) are the associated Laguerre 
polynomials, ri_o = [(2/2 - yi) - j(a;i - X2)]/(i\/2) and ro,i = [(yi - y2) - j(a;i - 
X2)\l {L\Fi). The gauge-dependent phase factor is 

X= (a;i-a:2)(yi+y2)/2L^ (29) 
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It is not translational invariant. The integrals over are 

^1 = ^ / 1? = ^ifo(|C|a«), (30) 



D2^^ dfc. = ^ifodCkn), (31) 

^3 - ^ / = ^ 7 + A{aA , (32) 

in which Z12 = zi — Z2 and 



a„ = ^Jl + rfn. (33) 

Further, Ac = Ti/{mc) is the Compton wavelength, C = (-^i ^ Z'i)l'^ci V = 

{2h(jJc/mc^y^^ , the cyclotron frequency is ujc = eB/m, and ^0(0;) is the Bessel 
function of the second kind. The quantity A{anC) is 



/ i^TT \ Z Z 

sin l^—j Iu{anO + 2^''(""C) - -j^JAf^nO 



(34) 



sin(7ri^)a„ 

where (i) is the Bessel function of the second kind in the standard notation and (t) 
is the Anger function |23j . 

Thus we expressed kernels K^(ri,r2) by combinations of S{ri — r2)l, (5(ri — r2)/3 
and G>o(t'i, T'2) operators, see Eq. (0). Operator 6o(t"i, ^2) of Eq. ([^ can be rewritten 
in a block form 

V T2(r-i,r2) -Ti{ri,r2) J 



where 



n— ^ ^ 



(36) 



71, n — 
•ii — y> 

Using Eqs. ([TT]) and ([20)1 we also can rewrite kernels A'^(ri,r2) in a block form 
where we denote (5 = 5{ri — r2)l. Then the kernels are 



A:*(ri,r2) 

/ (5(ri-r 

E 



2V^ 



n=0 



/ (5(ri-r2) (Da - iDi)A„_i^„_i - V^i:)2A„_i,„ 

(5(ri-r2) -V^L'2A„,„_i -(Dg + ii:)i)A„,„ 

(Dg + ii:>i)A„_i,„_i -V7^i:)2A„_i,„ (5(ri-r2) 

V -y^D2An,n-l {-D3+iDi)A„,n <5(ri - r2) 

Equation (jSH) together with Eqs. ([Ml) - (IMl), dSQ]) - (|32]) are the final expressions 
for the matrix elements of kernel operator. The kernels A'*(ri,r2) are given 
by combinations of two parts: translational-invariant gauge-independent parts and 
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the gauge-dependent phase factor e^^. Since x is not translational invariant, the 
kernels K^{ri,r2) are functions of independent variables ri and r2. 

The non-diagonal elements of matrix p9p are non-local in ri and r2 coordinates. 
In consequence, they change shapes of functions transformed with the use 
of K^{ri,r2), see Eq. pO)) . The transformed functions \^'^(ri) are more delocalized 
than the initial functions ^'(ri), since the non-diagonal elements of matrix (p9|) are 
characterized by the finite extent on the order of Ac or of magnetic length L. We 
describe this in the next section. 

The kernels K^{ri,r2) depend on two independent variables, i.e. on six 
coordinates. This case is more complicated than the field-free situation, where 
the transformation kernels depend on one independent variable R = ri — r2, see 
Refs. [niHn]- to slmpUfy our analysis, we set r2 = and denote ri = r. Using this 
simplification the matrix ([M)) can be written as 

( 5{r) f3(r)~zfi(r) ~f2(r) \ 

5{r) ^T2{r) -T3{r)-iT^{r 

f3(r)-Hzfi(r) -f2(r) 5{r) 

V -T2{r) -T^{r)+iTi{r) 5{r). J 

In the above equation instead of Di and A,„_„ integrals we defined another three 
functions 

ri(r) = \^^2 (41) 

n=0 

r2(r) ^ ^^^-"J ^2 E^o(«"ICI)LLi(p'), (42) 

^ n=l 



(40) 



r3(r) 



n— 



-+A(anC) 

z 



n— 

We introduced the notation — {x^ -f- y^) j (^l?) . In the second line of 
Eq. (j43l) we expressed the infinite sum over Laguerre polynomials by ^(p) function, 
see [Appendix B[ Functions ri(T) and r3(r) are defined similarly to ri(r), and r3(r) 
replacing an — )■ fln+i. The function r2(r) can be obtained from r2(T") replacing ri^o — > 
ro,i, see Eqs. ([27l) - ([28| . 

The summations over n in Eqs. (j4ip - (j43l) extend to infinity but in practice 
they are limited to finite numbers of terms by an exponential decrease of i4ro(a„|C|) 
and AiflnC) with n. However, for low magnetic fields there are large numbers of 
Laguerre polynomials involved in the summations so we truncate them for iiro(a„|^|) < 
10~^ and |y^(a„C)l < 10^^ or, alternatively, for n > 10^. This approximation allows 
us to perform calculations of Ti{r) for B > 10^ T, see below. For strong magnetic 
fields there is 77 — > 00, so the arguments a„C = {z / Xc) y^^+rfn of the Bessel or Anger 
functions in Eqs. (HT|) - P5|) are large for n ^ 0, and the terms including oqC do not 
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Figure 1. Calculated contour plot of ri(r), as given in Eq. II41I I. for the 
plane y = 0. Contour lines start for ri(7') = 0.05A^^ (outer lines) and end 
for ri(r) = lOAJ^. Coordinates x and z are in Ac units. For low magnetic fields 
the kernel is almost symmetric in x, y, z coordinates, while for large B the kernel 
in X and y directions is much more confined than in z direction. 



depend on the field. Then the summations over n reduce to single terms with n — 0. 
In this limit there is r2(r) ~ 0, while 



e 



-p^/2+ix 



ri(r)^ ttAL^ ^o(kl/Ac), (44) 



r3(r-) - + -^(^/Ac). (45) 



iS{p^) ie-P'/^+'>^ 
ttXcL'^z ttXcL"^ 

Thus both ri(r) and r3(r) have the extent L in the x — y directions and Ac in the z 
direction, so they are strongly anisotropic. For ri(r), r2('') and r3(T) functions there 
are no field- independent terms in the series (because oq — oi), so all these terms 
vanish in high magnetic fields and no widening occurs. 

In Figure[T]we show a contour plot of Fi (r) function for the y = plane calculated 
for various magnetic fields. For each of the four sections the contour lines start 
from Fi(r) = O.OSAj^ and they increase logarithmically to lOA^, ^. For all magnetic 
fields and all x and z the values of ri(r) are positive and they monotonically decrease 
to zero with growing x and z. For z = the expressions for Fi(r) diverge [see Eq. (PT|)]. 
but for X — and nonzero z they are large but finite. 

In Figure[T^ we show Fi (r) for a magnetic field B — 4.4x 10^ T, which corresponds 
to 1% of the Schwinger field of Bq = 4.4 x 10^ T. For such a field Fi(r) is, to a 
good approximation, spherically isotropic. The characteristic width of Fi(r) is on 
the order of Ac in all directions, as expected for the zero-field limit of the Dirac 
equation ^9j. Similar results are obtained for lower magnetic fields. At the Schwinger 
field (Figure [TJj) the term Fi(r) is still almost spherically isotropic but the widening 
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of ri(r) is larger than in the low-field limit. For still higher fields the terms ri(r) 
become anisotropic: their extent in the z direction changes slowly but they significantly 
shrink in the x — y directions. The behavior of ^2{f) and r3(r) functions is similar to 
that described for Fi (r), however F2(r) vanishes in the limits of i3 — >■ oo or i? — 0. For 
magnetic fields lower than the Schwinger field the behavior of Fj(r) terms is similar 
to the analogous Fj(r) terms. For higher fields all these terms go to zero. 



3. Second moments of transformation kernels 



Now we estimate quantitatively the spatial extent of the elements of kernel 
matrices K^{r,0) given in Eq. pO)) by calculating second normalized moments Wj 
of Fj(r) and Tj{r) functions (j=l,2,3). Here we follow Rose who used second 
moments to characterize non-locality of the transformation kernel for the Foldy- 
Wouthyusen transformation in absence of external fields. The second normalized 
moments are defined as 

/ r'^rj{r)d^r 
/r,(r)d3r 

The second normalized moments Wj of ^j{r) functions are defined similarly. The 
above integrals over F2(r), r3(r), F2(r) and F3(r) functions vanish due to symmetry. 
Below we analyze in more detail the second moments Wi and Wi of the ri(r) 
and ri(r) functions. Introducing auxiliary quantities 

Mt^^ f z''Ko{\z\aJ\c)dz, (47) 



where p = p/{-j2L), we obtain 



(48) 



2 kAP 

2 /„3 



Direct calculations give A^g — '^/o-n and A^f = A^/a„, where a„ are defined in 
Eq. ([33| . On setting in Eq. pS]): v = Q and t = p^, making the substitution xy/ (2L^) = 
sin(20)/(4L^) and integrating over d(j), we obtain 

/•oo 

Xg= / e-*/Vo(t/2)L°(t)dt - \/2P„(0), (50) 
where P„(a;) are the Legendre polynomials. Using the above results, we have 

f:A^gA^s=v2f:^, (51) 



n— 



Y^MlM^,^V2XlY.-^. (52) 
To calculate proceed in a similar way and obtain 



= 2L2 / e-*/Vo(i/2)tL°(i)dt. (53) 
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In order to use the result ([50)) we must express the product in terms of L°„(i) 

polynomials. Using the identities: (n + l)L^_^_^{t) - {2n + 1 - t)L° (t) + nL^^_^{t) = 
and L^{t)t = L^{t) - L^(t) we have 



or2 /-oo nr2 poc 

^ = / e^''^Jo{t/2)hl{t)tdt + / e-'/^Mt/2)hl{t)tdt 

f— ; an Jo o-o Jo 



n=0 n—1 



2i^V2<' Vp„(0) {---^) + VnP„(0) _ ^ _ ^ ) ^ , (54) 



I \0.n Q-n+l / — „ \CI71 0,n+l Q-n-l 

where we put l/a_i = and used the identities: P2n+i(0) — and nP„(0) = for n = 
0, 1. The sums over n can be calculated numerically either noting that for m — 2n there 

is Pm(0) — (—1/4)™ ( 1 or, alternatively, by a direct summation with the use of 



the generating functions for the Legendre polynomials, see [Appendix A Calculating 



the second moment of f i(r) one should substitute a„ a„+i in equations ([5T[) . ([52[) 
and ([54[). 

In the limit of high fields there is r/ ^ 1 and l/a„ ~ Sn^, so the second 
moment of ^i{r) in the x — y directions is Wi ~ 2L^ cx B~^, while the second 
moment of ri(r) in the z direction is Wf ~ A^, see Eq. ([44[) . For Wi there are 
no field-independent terms in the summations in Eqs. ([FT[) . (15^ and ([M|. and in 
high fields there is VV^ ~ Wf ~ 0. Thus no widening occurs. For low magnetic 
fields, when 77 ^ 1, we may approximate a~P ~ 1 ~ prfn. Since X]nPn(0) — 
(see [Appendix A I we have from Eq. ([5T[) that lims^o X^^^o -^o-^o = 1- This gives. 



upon using Eqs. ([49|), ([52|) and ([54[) . 

00 

lim ~ \/2A^ ^ P„(0)(1 - 3772^/2) ~ A^, (55) 

ra=0 

lim Wf ~ 2L2v^ <^ V P„(0)^ \ = 2\l, (56) 



The same results are obtained for Wf and moments. Thus the second moments 
of ri(r) and V\(r) functions in the low-field limit are 

lim Wi = lim Wi = 3A2. (57) 

This result reproduces our value obtained previously for MOT at _B = 0, see Ref. |10) . 

In Figure [2) we show Wi (thick line), dashed line) and Wf (dash-dotted line), 
as defined in Eqs. (HHl) . ([Fl]) . ([5^ and ([M[) . versus magnetic field B. For low magnetic 
fields WC tends to 2A^ and remains nearly constant up to the Schwinger field, where 
a flat maximum occurs. For larger magnetic fields starts to decrease and finally 
decreases as . The second moment remains nearly constant, Wf ~ A^, within 
the whole range of magnetic fields except in the vicinity of the Schwinger field. The 
total second moment Wi = Wf -t-Wf of the ri(r) function, exhibits for an increasing 
magnetic field a transition from the 3D value Wi = 3A^ to the ID value Wi — A^. 
For low magnetic fields and Wf moments are close to and . However, for 
magnetic fields higher than the Schwinger field both of them tend to zero exponentially 
with B. It should be reminded that the kernels given by Eq. P0[) are composed of 
various elements. All of them are nonlocal (with the exception of the delta functions 
on the diagonal), but only ri(r) and ri(r) have non- vanishing second moments. It is 
the second moment of ri(r) that is shown in Figure [2) 
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Figure 2. Calculated normalized second moment of ri(7'), as given in 
Eq. I|49| l. versus magnetic field B. Solid line: total second moment Wi, dashed 
line: Wj , dot-dashed line: Wj . Second moment Wi exhibits a transition from 3D 
case Wi ~ 3A|? at low fields to ID case Wi ~ at high fields. 



For a non-relativistic electron in a magnetic field there is no widening of the 
transformed wave function. We retrieve this limit by setting mc^ — > oo, which 
gives: — > 0, — > and a,i = 1 for all n. Then the sums in Eq. vanish, 
the summation in Eq. (jSip gives unity, and the summation in Eq. (j52p gives — > 0. 
Therefore = 0, VVf = 0, VVf = and no widening occurs. 



4. Transformed functions 



In this section we consider properties of transformed functions using as an example 
a Gaussian function of finite width. We take the initial function 'i'{r) in the 
form: ^'(r) = /(r)(0, 1, 0, 0)^, in which 

/M = -Vid^' (58) 

and d characterizes function's width. The transformed function is: ^'{r) = V'i>{r) = 
/(r)/v/2(0, 1,0, i)^- Then we have l^-'*) = which gives [see Eq. ^] 



V2 



/ 



± 



2 ^ 

n=0 




firi 


-iy^D2{zi - Z2)A„_l,„(pi, P2) 
{£'1(21 - Z2) - iDsizi - Z2)}An,„{pi,P2 

-y/nD2{zi - Z2)A„_i,„(pi, P2) 

\ {iDi{zi- Z2) - D3{zi- Z2)}An,„{pi,P2) J 



/(r2)d'r2,(59) 
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where A,„ „ and Di are defined in Eqs. - and ([50]) - ([5^ . As seen from 
Eq. ([5^ . positive or negative energy functions ^'^(ri) consist of two parts: the 
function ^''(ri) and the nonlocal parts arising from the integration of ^''(ri) over 
appropriate elements of r2). The summation over n and the threefold integrals 

over r2 can be performed numerically. 

Now we turn to the high-field approximation, in which L ^ A, 
the summation over n in Eq. 
transformed function is 



59|) reduces to the single term with n 



In this case 
= and the 



*'±(r) 



/2C 



/ 








/ 






f{r) 


± - 






Fz{z)Fp{p) 







2 







V 








V 


iF,{z)Fp{p) j 



where C is a normalization constant and 

/2d /"^ (1 - acfc^)e-'=''*'/2+''=^ 



27r3/4 



The presence of Gaussian term ensures the convergence of integration. Further 

1 



(60) 



(61) 



Fp{p) 



2V27rd 



P2 



V(F+2(Pl2 + 2l4 



-a^x —a-yV +«o xy 



where 



2(d4 + 2d2L2 ^ 



2L4)' 



6^ = 



2L2(d4_ 



2d^L^ + 2L4) ^ 



d4 + 2^2^2 + 2X4' 



For narrow initial functions: d L, there is 



(62) 

(63) 
(64) 
(65) 

1/(4L2) and b"^ ~ l/(2£2), 



so that Fp{p) has a spatial extent on the order of L/2, similar to the case analyzed 
before, see Eqs. (|311) and P5|) . In the opposite limit of wide initial functions: d ^ L, 
we have ~ &2 ^ 1/(2^^) and Oy ~ 1/{2L'^). The differences between Ox and Cj, 
result from the asymmetry of the Landau gauge. To obtain the same results in another 
gauge, one would have to modify the initial functions accordingly, see Refs. [26l [27] . 
To demonstrate a widening of the functions 'i''^{r) we calculate their variances 



= j 'i!'^{r)^rH'^{r)d^r j J \^'^{r)\^d^r , 



(66) 



which are different from the second moments calculated above, see Eq. (|46p. The 
integrals in Eq. can be calculated analytically but we do not quote the long 
resulting expressions. 

The plots of for various magnetic fields and initial function widths are given 
in Figure [3l The solid lines represent calculated variances of the transformed function. 
The variances Zc = Sd'^ /2 of the initial functions are indicated on the right ordinate. 
For the initial function of the width d the variances of the transformed functions are 
larger than Zq, so the transformed functions are wider than the initial one. For the 
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Figure 3. Calculated variances of transformed Gaussian function ^'^(r), as 
given in Eq. Il60|l . versus magnetic field B/Bq in the high-field approximation. 
Solid lines: variances corresponding to the initial functions having different 
widths d. Variances of the initial Gaussian functions for d = Xc and d = 0.5Ac 
are indicated. 

initial functions having d ~ Xc the variances reach their limit Zq for magnetic 
fields on the order of = 10^^ T. For smaller d, the variances of transformed functions 
are much larger than Zq and they reach the limit Zq for very large fields. All results 
presented in Figure [3] are valid for magnetic fields larger than the Schwinger field, for 
which L ^ Ac- The asymmetry between positive and negative variances Z^ shown 
in Figure [3] arises from the asymmetry in the initial wave function having only one 
nonzero component. Selecting other nonzero component or more than one nonzero 
components, one would obtain different numerical values of Z^, but their general 
properties would be similar to those presented in Figure [3l 

5. Tv^ro-dimensional case 

Recently, Lamata et al. j21| proposed new topological effects for the two-dimensional 
Dirac equation in a magnetic field. They also indicated a possibility of simulation 
and reconstruction of the Wigner function for the lowest Landau levels with the use 
of trapped ions. In this connection it is of interest to calculate the non-locality of the 
transformation kernel in two dimensions (2D). The matrix elements of 2D kernels can 
be obtained from the three dimensional results given above by setting e'*^^^ — )■ 6{kz) 
in all integrals including kz and by a proper adjustment of normalization constants. 
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Figure 4. Gaugc-indcpcndent part of kernel's element Gi(p) in two dimensions, 
as given in Eq. JTO}, calculated for four values of magnetic field. For low fields the 
kernel non-locality is on the order of Ac and it does not depend on the field. For 
higher fields the non-locality slightly increases and for very large fields it decreases 
again. 



The 2D counterpart of Eq. gOl) is 

^/3± i 



/ S{p) 

6{p) 

+iGi{p) -G2{p) 

\ -G2{p) +iGiip) 



-iG^{p) 

~G2(p) 

5{P) 




-iG,{p) 


<5(P). / 



(67) 



In the above expression we defined 
Gi(p) 



G2{P) = 



L2 



E 



(69) 



using the same notation as in Eq. (|40p . Functions Gi{p) and G2{p) are defined 



similarly to Gi{p) and G2{p) replacing 



rl^{n + 1), see Eqs. ^ 



Note that in the 2D case there is no counterpart of the r3(r) term. 
Functions Gi{p), G2(p), Gi{p) and G2{p) are products of gauge-independent parts 
and the gauge- dependent factor e'^. The summations over n in Eqs. (|68t - (|69| 
are slowly convergent. However, these sums can be transformed to fast converging 
integrals that can be calculated numerically. For Gi{p) one has 
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where t = exp(— 77^(7^). Details of the derivation are given in [Appendix A[ The 
remaining sums can be transformed into integrals in a similar way. 

In Figure |3] we plot the gauge- independent part of Gi{p) for four values of 
magnetic field. For low fields the spatial extent of the kernel is on the order of Ac 
and it does not depend on the field. For the Schwinger field (dash-dotted line) the 
non- locality of Gi{p) differs only slightly from its low-field value (solid line). For 
larger fields the non-locality of Gi{p) increases (dashed line) and for very large fields 
it decreases again (dotted line). The non- locality of Gi{p) is on the order of the 
magnetic length L which in this regime is much smaller than Ac. For large fields in 
three dimensions the non-locality of the kernels tends to zero in the x and y directions, 
but it remains finite in the z direction. For large fields in 2D the non-locality of the 
kernels goes to zero. Therefore, for large fields in 2D there is no widening of the 
transformed functions. This conclusion may be important for simulations of the Dirac 
equation by trapped ions 12 Ij . where arbitrary large values of a simulated magnetic 
field may be achieved. 

6. Discussion 

As follows from our considerations and, in particular, from the above figures, the 
magnetic field effects for relativistic electrons in a vacuum become appreciable at 
gigantic field intensities, comparable to the Schwinger field Bq — 4.4 x 10^ T. 
However, one can overcome this problem in two ways. The first is provided by 
narrow-gap semiconductors, whose energy band structure can be described by Dirac- 
type Hamiltonians with considerably more favorable parameters, see Ref. [151 • For 
example, in InSb the energy gap, corresponding to 2mc^ in a vacuum, is Eg ~ 0.23 
eV and the effective mass m*, corresponding to rest electron mass m in a vacuum, 
is m* ~ 0.014rn. In consequence, an effective Schwinger field in InSb, given by the 
condition heBQ/m* — Eg/2, is Bq ~ 14 T. Such magnetic fields are easily available 
in laboratories, so the field effects described above can be readily investigated in 
narrow-gap semiconductors. One should, however, bear in mind that in experiments 
with crystalline solids one deals with additional effects, most notably with electron 
scattering leading to broadening of Landau levels, which were not taken into account 
in our considerations. 

The second way to lower effective magnetic fields is to use simulations of the 
Dirac equation employing trapped ions or cold atoms interacting with laser radiation, 
see e.g. Refs. |T71 [THl [H] and review [7]. Such simulations may also include the 
presence of an external magnetic fields [2j [21] . In this type of experiments an effective 
simulated field can be tailored to suit observation requirements. For example, in a 
simulation of the Dirac equation proposed in Ref. [2] , the effective simulated magnetic 
field corresponds to the ratio (fieB / ra) l(2mc^) = 16.65 for realistic trap parameters. 
However, in simulation- type experiments there is no need to apply real magnetic fields 
apart from the fields necessary to create ion traps. 

As to our results, the general rule is that, similarly to the non-locality of 
energy separating transformations for i? = 0, the characteristic quantity is the 
Compton wavelength Ac = h/mc, see Ref. [lO]. In the presence of a magnetic 
field, typical anisotropy appears between the z direction along the magnetic field 
and the transverse x — y plane. This is seen in the broadening of the transformation 
kernel illustrated in Figure [1] The anisotropy appears also in the second moments of 
functional transformation kernels shown in Figure [21 This result can be considered to 
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be a direct extension of that of Rose [5], who characterized quantitatively the non- 
locahty of the Foldy-Wouthuysen transformation for B = 0. The low field limit Wi = 
3A^ agrees with our previous calculation for MOT at B = 0. The corresponding result 
for FWT is Wi = (3/4) and one can say that the FW transformation is somewhat 
more "compact" than the MO transformation. Interestingly, at high fields Wi is 
suppressed, which results in Wi — . This means that at high B the quantization of 
the electron motion in the x—y plane gives in consequence an effective one-dimensional 
motion along the magnetic field. The parallel second moment Wf depends somewhat 
on B because in the relativistic mechanics there is no exact separation between parallel 
and transverse directions of magnetic motion. Finally, Figure [3] indicates two effects. 
First, the broadening of transformed functions is stronger for narrow initial widths. 
Second, with an increasing magnetic field the broadening of transformed functions 
diminishes, reaching asymptotically the non-broadened values at very high fields. One 
should bear in mind that the results shown in Figure [3] were obtained in the high field 
approximation, i.e. B / Bq > 1. 

The results obtained in this paper refer to a non-locality of the energy-separating 
transformation for the Dirac Hamiltonian in the presence of a uniform magnetic 
field. However, under two assumptions given below it is possible to generalize our 
approach to the case of an inhomogeneous magnetic field. First, there should exist 
an exact transformation separating the positive and negative energy states. Second, 
it is practical to have analytical solutions of the Dirac equation for electrons in an 
inhomogeneous magnetic field that allow one to find an expansion of the kernel into 
eigenvectors in a way analogous to that given in Eq. ([T2l) . More specifically, the 
existence of transformation separating positive and negative energy states can be 
simplified to an assumption that the operator Ti ^ [see Eq. Q], describing the square 
of the Dirac Hamiltonian in an inhomogeneous magnetic field, is diagonal. This 
criterion can be satisfied for example by a vector potential having the form A = 
[Ax{x,y), Ay{x,y),0] with arbitrary functions Ax{x,y) and Ay{x,y). Then the 
operator separating the positive and negative energies is analogous to U given in 
Eq. ([5]) but with a suitably changed vector potential. Physically, the existence of such 
an operator is limited to situations in which there is no creation of electron- hole pairs. 
Therefore, such an operator does not exist for Dirac particle in an electric field or in 
some configurations of the magnetic field that mix positive and negative energy states. 
As to the second requirement for the applicability of the present approach, there are 
several inhomogeneous magnetic fields known in the literature for which analytical 
solutions exist. As an example, in three dimensions the Dirac equation has solutions 
for vector potentials A = [(B/a) tanh(ay), 0, 0] or A = [{B/a)[l — exp{ay)],0,0], 
where is a characteristic length of the potential |28j . In two dimensions there 
exist also solutions for magnetic field B = /3 sec^ {ay) . The above vector potentials 
have the form A — [Ax{x, y),Ay{x, y),0], so that the operator H ^ is diagonal and one 
can calculate the non-locality of the transformation kernel in a way analogous to that 
presented above. We expect the results to be similar to those presented in our work 
with properly modified magnetic length. 

Our calculations of the non-locality generated by FWT and MOT at S = 0, see 
Ref. [lOj . indicate that different "separating" transformations have similar properties 
and small differences between them concern mostly the extent of non-localities. 
In consequence, we believe that the properties of MOT investigated in this work 
characterize other separating transformations for the Dirac electron in the presence of 
a magnetic field. 
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7. Summary 

We investigate non-locality of transformations separating positive and negative energy 
states of the Dirac electrons in the presence of a magnetic field, using as an example the 
Moss-Okninski transformation. Functional kernels generated by the transformation 
are described and their second moments in the coordinate space are calculated to 
determine quantitatively the nonlocal features. The latter are characterized by the 
Compton wavelength Ac = h/mc. The behavior of kernels is different along the 
magnetic field and in the plane transverse to it. Dependence of non- locality on the 
intensity of magnetic field is studied. It is demonstrated that the transformed functions 
are broadened. Finally, it is indicated that a practical way to study Dirac electrons in 
the presence of a magnetic field is to simulate their properties using either narrow-gap 
semiconductors or trapped ions interacting with the laser radiation. 



Appendix A. 

Calculating integrals Am,n in Eqs- ((26|) - p8| we use the identity 

/oo 
H„,((7 + a)H„((j + b)e-'''dq = 2"V^m!6"'L;^' (-2a&), (A.l) 
-oo 

where n' = n — m and m < n. 

A convenient way to calculate sums in Eq. ()54|) is to use the generating 
functions for the Legendre polynomials: X^i^o — ~ ^te -I- 

and X^i^o ~ ~ ~ 2te -1- t^)'^/^. In order to apply the generating 

functions we first make the transformation l/a„ = 2/y/TT exp[—q^ — q^ri^n]dq, next 
we set t = e~'^ , then change the order of summation and integration. As a result 
one obtains 

^ Pn(0) - 2 e-o'dq 

nP„(0) r e'^'-^^'^'dq 

For vanishing 77 (weak magnetic fields) the above integrals tend to l/\/2 and — 1/\/8, 
respectively, while for 77 1 (strong magnetic fields) they tend to ±1. 

The same method can be used to calculate a sum of the type Xn-Pn(0)/(1 + 
Ty^n)'^/^, employing the identity 1/(1 -I- r/^n)'^/^ _ —2[d/db)\b-\-Tfn]~'^/'^ taken at = 1. 
Proceeding the same way as for Eq. \A.2\ we get 

P„(0) ^ J_ r g^e-^^dg 



which applies to Eq. (|52l). 

Consider the sum in Eq. (|68| for the 2D case. Using the identity: 1 / -^/l + rfn = 
(1/0F) e-«'-™''«'dg we obtain 

s = f; = 4 i: /" ^sf)'-'' U-''")' 
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Changing the order of summation and integration, and using the identity 

flKi^r = ^^3^exp (A.6) 

with t — exp(— 77^(7^) one obtains Eq. ([70)1 . The presence of a Gaussian term in 
Eq. (|A.5I) ensures fast convergence of the integral for large q. For smaU q the 
factor exp[te/(t — 1)] tends to zero as exp(— l/q^) and the integrand in Eq. (155)) 
remains finite for all values of q. 



Appendix B. 

We prove a sum rule for the A„^„ integrals. Let pi and p2 be position vectors in 
the two-dimensional space and {|q)} the complete set of eigenstates of a hermitian 
operator. Then 

^(Pl - P2) = (P1|P2> = 5](pi|q>(q|p2> = 5] Vq(Pl)V'q(P2)^ (B.l) 

q q 

where V'q(p) = (pIi)- Taking the eigenstates |q) = \n,kx) of the Hamiltonian for a 
non-relativistic electron in a magnetic field we obtain 



/ on ^ rr H„(a)H„(e2)dfc.-<5(pi-p2), (B.2) 

where p = {x,y). Performing in Eq. (jB.2p the integration over kz and using the 
notation from Eq. (|25|) we have 



E 

ra=0 



n.n 



5{pi - P2). (B.3) 



Identity (jB.3[) was used in Eq. for the calculation of infinite sums over the Laguerre 
polynomials appearing in the l/z term. 
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Abstract. We investigate a non-locality of Moss-Okninski transforma- 
tion (MOT) used to separate positive and negative energy states in the 3+1 Dirac 
equation for relativistic electrons in the presence of a magnetic field. Properties of 
functional kernels generated by the MOT are analyzed and kernel non-localities 
are characterized by calculating their second moments parallel and perpendicu- 
lar to the magnetic field. Transformed functions are described and investigated 
by computing their variances. It is shown that the non-locality of the energy- 
separating transformation in the direction parallel to the magnetic field is char- 
acterized by the Compton wavelength Ac = h/mc. In the plane transverse to 
magnetic field the non-locality depends both on magnetic radius L = (h/eB)^^^ 
and Ac. The non-locality of MO transformation for the 2+1 Dirac equation is also 
considered. 
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1. Introduction 

The Dirac equation, in spite of its numerous applications and fundamental significance 
in relativity and quantum mechanics, is far from being completely understood. This 
is especially true for the Dirac electrons in the presence of external fields. For 
example, it has been known for many years that the Dirac equation can be solved 
exactly in case of an external uniform magnetic field [1], but it is only recently that 
an explicit electron dynamics of Dirac electrons in the presence of a magnetic field 
was worked out [2]. This problem is related to a somewhat mysterious phenomenon 
known in the literature under the name of Zitterbewegung. The phenomenon was 
conceived in 1930 by Schrodinger who observed that, according to the Dirac equation, 
operators corresponding to the electron velocity do not commute with the Dirac 
Hamiltonian, so the velocity depends on time also in absence of external fields [3]. 
It is understood at present that this remarkable phenomenon is due to an interference 
of states corresponding to positive and negative electron energies [31 [6]. In fact, it 
is known by now that, whenever one deals with a spectrum of positive and negative 
energy branches, the Zitterbewegung will appear [7|. 

There exist ways to avoid the duality of positive and negative energies. The best 
known of them is the Foldy-Wouthuysen transformation (FWT) which, in case of no 
external fields, allows one to break exactly the 4x4 Dirac equation into two 2x2 
equations for positive or negative energies [S]. In their pioneering paper, FW remarked 
that a kernel transforming functions from the original to the FW representation is 
characterized by a non-locality in coordinate space on the order of the Compton 
wavelength Ac = h/mc. Rose in his book W put this statement on a quantitative 
basis by showing that the second moment of the functional kernel is equal to (3/4)A^. 
This result was not followed by other investigations and it is by now not well known. 
The problem was recently taken up by the present authors who demonstrated that also 
other transformed quantities are characterized by non-localities on the order of Ac [10] . 

Case 11 showed that also the Dirac equation including the presence of an external 
uniform magnetic field B can be exactly separated into two equations for positive and 
negative electron energies. This form allows one to find easily the energy eigenvalues 
in any desired semi-relativistic approximation (see our Eq. ([7]) and Ref. [51]). It is 
possible to separate the positive and negative electron energies for i? 7^ because the 
presence of a uniform magnetic field does not induce electron transitions between the 
two parts of the spectrum. The case of magnetic field is of interest for several reasons. 
First, it introduces an adjustable external parameter B affecting all energies and other 
electron characteristics. Second, it quantizes the orbital motion in two degrees of 
freedom. Third, it quantizes the spin energies. However, to our knowledge properties 
of the separating transformations for this important case were not investigated until 
now. 

In the form proposed by Case the transformation operator is C/ = e^^^'^/'^ with p2 
being the standard Dirac matrix, and the phase given hy (f) — tan~^(o'7r)/mc, where & 
are the Pauli matrices and tt is the generalized momentum operator including the 
vector potential of the magnetic field. In this form it is difficult to analyze the 
properties of U . A method to get a more explicit form of U was described by de 
Vries [12]. However, it was shown that the "separating" transformations, both in 
absence of fields and for B 7^ 0, are not unique. In other words, it is possible to 
separate the Dirac equation into two equations for positive or negative energies, 
both in absence of fields and in the presence of a magnetic field, using different 
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unitary transformations [T^ [T51 UM [T5] . In our description involving Dirac 
electrons in a magnetic field we use a two-step transformation proposed by Moss and 
Okninski (MO) [H]. An advantage of this choice is that the MOT is given explicitly, 
it is convenient for kernel calculations and we can compare our final results in the 
limit i? — > with our previous description for _B = 0, see |10| . In our treatment we 
do not concentrate so much on the separation of energies but rather on the properties 
of transformation operators and the transformed functions. Naturally, we emphasize 
the effects of an external magnetic field. 

Our subject is of relevance for at least three reasons. First, with the comeback of 
interest in narrow-gap semiconductors and, in particular, the appearance of zero-gap 
graphene [16], there has been a real surge of works concerned with relativistic-type 
wave equations. Second, it is now possible to simulate the Dirac equation with the 
trapped ions or cold atoms interacting with laser radiation, where one can tailor 
"user friendly" values of the basic parameters mc^ and c [H [T71 [TSJ [H]. In fact, a 
proof-of-principle experiment simulating the 1-1-1 Dirac equation and the resulting 
Zitterbewegung was carried out by Gerritsma et al. [20]. Third, one can simulate 
equally well the Dirac equation in the presence of a magnetic field [21 [21] . In such 
simulations it is possible to tailor a seeming intensity of the field making its effects on 
relativistic electrons relevant in terrestrial conditions. Putting it in more quantitative 
terms, one can achieve experimentally the so-called Schwinger field, defined by the 
equality heBo/m = mc^, not only in the vicinity of neutron stars but in standard 
laboratories, see our discussion below. 

Our paper is organized in the following way. In Sections II and III we introduce the 
Moss-Okninski transformation, describe its functional kernel and calculate its zeroth 
and second moments. In Section IV we transform Gaussian functions and calculate 
their normalized variances. In Section V we analyze kernels in two dimensions. In 
Section VI we discuss our results. The paper is concluded by a summary. In two 
Appendices we give some mathematical details of our derivations. 



2. Transformation Kernel and its non-locality 



In this section we define the MO transformation and its functional kernel. Then we 
analyze kernel's properties and, in particular, its non-locality. The Dirac Hamiltonian 
for a relativistic electron in a magnetic field is 

"H ~ COixTTx + CCtyTTy + CCtzT^z + ^mC^ , (1) 

where tt = p + 2^ is the generalized momentum, |e| is the electron charge, A is the 
vector potential of a magnetic field B, and di and /3 are the Dirac matrices in the 
standard notation. Following Moss and Okninski [15) we first introduce a unitary 
operator 

(2) 



in which 



/ 







' \ 











—i 


i 











\ 


i 





J 



(3) 
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This operator transforms the initial Hamihonian into a form with vanishing diagonal 
terms 

= caxTTx + cayTTy + ca^TTz + Smc . (4) 
The second step transforms "H to a strictly diagonal form with the use of an operator U 

?7 = (/3 + e) , (5) 

where 

e = „^ , . (6) 

(H'2)l/2 

After the complete transformation the Hamiltonian has a diagonal form (cf. Refs.|lll 

UnU^ = Pim^c^ + c^TT^ - ec^aSB)!/^, (7) 
where S = ( ^ ) is the 4x4 spin operator. Since tt is an operator, the square 



X cr ^ 

root is understood as an expansion in powers of tt. In the following we concentrate 
on the properties of the transformation and transformed functions. 

Let l^") be a function in the old space and \^') — V]"^) in the transformed space. 
We introduce standard projecting operators for positive and negative energy states 

F± = i(i±e). (8) 

If the projection operators are applied to the transformed functions l^*'), they split 
them into states corresponding to positive or negative energies: |^''^) — P^\^'). 
Next l^*'^) are transformed using the operator t/, i.e. |^'^) — UY^'^). We introduce 
operators IJ^ defined as = JJP^ with the property 1^'=*=) = IJP^\^'). We have 

c>± = f /3 + e) f i ± e 



2^2 

Let us consider transformed functions \f'^(ri) = (ri|?7*|5''). Inserting the identity 
operator 1 = / \r2){r2\<P'r2 one has 

k^{rur2W{r2)d?r2, (10) 

where we defined the functional transformation kernels 

k^{ri,r2) = {ri\U^\r2). (11) 

In the following we investigate the properties of K^{ri,r2)- To calcu- 
late K'^{ri,r2) we insert twice the unity operator 1 = |m)(m|. The states |m) 
are eigenstates of the Hamiltonian Ti, given in Eq. ([4]) , and the summation is carried 
out over all quantum numbers characterizing |m). Thus 

i^±(ri,r2) - Yl (n|ni)(m|;7±|m')(m'|r2). (12) 
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The operator consists of operators: i, (5 and Q, see Eq. The matrix elements 
of /3 and 1 are diagonal: /3m, m' — PSm.m' and l(5m,m', respectively. The non-diagonal 
contribution to f{^{ri,r2) arises from 9 

e(ri,r2) - (ri|m)(m|e|m')(m'|r2). (13) 

m,m' 

The eigenstates |m) of the Hamiltonian H are related to the eigenstates |n) of the 
initial Dirac Hamiltonian H by a unitary transformation 

|m) = V\n), (14) 

where V is defined in Eq. ([2]). We want to calculate the elements of 0(T'i,r2) matrix 
and for this purpose we need to know explicit forms of the eigenstates |m). 

We take the vector potential of a magnetic field in the Landau gauge A = 
{—By, 0, 0). Then the eigenstates |n) depend on five quantum numbers: n, kx, kz, e, s, 
where n is the harmonic oscillator number, kx and are wave vector components, e = 
±1 labels positive and negative energy branches, and s — ±1 is the spin 
index corresponding to the spins ±1/2, respectively. In the Johnson-Lippman 
representation [22] the eigenfunctions ^ni^) = (''In) of the Hamiltonian H are 

f si{eEn^k, + mc^) (t>n-i{r) \ 
S2{eEn^k, + mc^) (/>ri(r) 

(siCPz - S2huJn) (/>«-i(t-) 
\ -{sihuJn+ S2CPz) (/>n(r) J 

where si = (1 + s)/2 and S2 = (1 — s)/2 select the states s = ±1, respectively. The 
frequency is aj„ = ujy/n with uj = \/2c/L, the energy is 



nepz 



(15) 



En,k^ = V(mc2)2 + (nw„)2 + (UkzC)^, (16) 
and it does not depend explicitly on s. The norm in Eq. ((T5|) is Nnek^ — i'^^n k + 
2emc^£'ri,fe^)~^/^, and the (/>„ functions are 

Mr)^ 2^7ZC ^-^^^^ ' 



where H„(^) are the Hermite polynomials, C„ = •\/2"n!^/7r, the magnetic radius 
is L = \fhJeB and ^ = yjL — kxL. 

Since V defined above is an unitary matrix, there is -Em — E^ — eEn^k^- The 
states |m) are characterized by the same five quantum numbers describing |n), so there 
is also (5m, m' = ^n,n'- For two eigcnstatcs |m) and |m') of the Hamiltonian T-L we have: 
{m\n'\m') = Em(5m,m', {m\H'^\m') = EfJ^^n,', and 

E 

(m|e|m') = (5m, m' = qeSnAi', (18) 

\/ E 



where q — 1 or q = —1. Selecting q = 1 we obtain from Eqs. (fT5)) . and (HI 

(nieira) =^e(ri|m)(m|r2) =^e(ri|\/n)(nf|r2). (19) 

m n 

Selecting = — 1 in Eq. would not have changed the results because we perform 
a summation over e — ±1. Since 1^ is a number matrix, there is (r|yn) = V{r\n), and 
one can rewrite (t'i|6|t2) in the form 

{ri\e\r2) =Veoiri,r2)V\ (20) 
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where we introduced an operator Qq by the relation 

(ri|eo|r2)=^6(ri|n)(nh). (21) 

n 

The advantage of 6o('"i, '"2): as compared to 6(ri,r2), is that 6o('"i,'''2) can be 
directly expressed in terms of the eigenstates |n) given by Eq. (ITS]) . Writing in Eq. (PT|) 
all quantum numbers we have 

/■oo /"OO 

(ri|eo|r2) = / dkj dfc. ^ e$n(n)$l;(r-2), (22) 

where $n(^) are given in Eq. (|15p . The integrals over and as well as the 
summations over e and s may be done analytically, the summation over n is to be 
performed numerically. However, in the limit of strong magnetic fields analytical 
results may be obtained. The summations and integrations indicated in Eq. ()22|) are 
carried out in the following order. First, we sum over the spin variable s noting 
that S1S2 = and J2s^i ~ 12s^2 — 1- I^i the summation over e we use the 
property: J2e=±i e = 0. As a result, one obtains 



poo I* 

eo(ri,r2) = ^ / dkJ 



dk. 



,(23) 



eo(ri,r2) = £ 



/ 'yi(f>n'{ri)<Pl,{r2) 73'/>n'(ri)<?!);,(r2) -720n' (r-i)*/),* (r2) \ 

7i(/)„(ri)(/),*(r2) ~720n(r-i)C'('^2) -730n(ri)0;(r2) 

73(f>n'{ri)<j)*^,{r2) -720n'(n)0;(r-2) -Ji(l)n'{ri)(f>*^,{r2) 

V -I2<l>n{ri)<l>n'{r2) -73 (n )(/),* (r2) -7i0n(ri)C(r2) / 

where 71 = mc^ /{2ttE), 72 = ?icij„/(27riJ), 73 = cpz/i^nE) and = n — 1. The 
matrix elements of 0o(t'i, r2) are products of integrals over k^ and k^. There appear 
three different integrals over k^: Di, D2 and D3 and four integrals over k^ denoted 
as Am,n- We write explicitly 

^ i:>lA„_i,„_l D^An-l^n-l -y/nD2An-l.n \ 

n,n — 

£'3A„-1,„-1 ~^D2An-\,n ~DxAn-\^n-\ 

\ -y/nD2An^n-\ -D-iAn^n -DiA„_„ / 

The integrals over k^ are 

Am.n= / e-M"^ 7^ H„,(a)H„(g2)e'^/'^^'+^^^dfc., (25) 
where = yj / L — k^L with j = 1,2. They are calculated analytically 

A„-,,„-, = ;^L^^.(p^)e-'^V2+.x^ (26) 

A„_i,„ =-^L^,(p-^)e-'^^/^+^^ (27) 

A„,„_i =^L^,(p-^)e-P^/2+»^ (28) 

where j = 0, 1, p = [(a;i — 2:2)^ + (yi — ?/2)^]/(2i^), L^(p) are the associated Laguerre 
polynomials, ri_o = [(2/2 - yi) - j(a;i - X2)]/(i\/2) and ro,i = [(yi - y2) - j(a;i - 
X2)\l {L\Fi). The gauge-dependent phase factor is 

X= (a;i-a:2)(yi+y2)/2L^ (29) 
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It is not translational invariant. The integrals over are 

^1 = ^ / 1? = ^ifo(|C|a«), (30) 



D2^^ dfc. = ^ifodCkn), (31) 

^3 - ^ / = ^ 7 + A{aA , (32) 

in which Z12 = zi — Z2 and 



a„ = ^Jl + rfn. (33) 

Further, Ac = Ti/{mc) is the Compton wavelength, C = (-^i ^ Z'i)l'^ci V = 

{2h(jJc/mc^y^^ , the cyclotron frequency is ujc = eB/m, and ^0(0;) is the Bessel 
function of the second kind. The quantity A{anC) is 



/ i^TT \ Z Z 

sin l^—j Iu{anO + 2^''(""C) - -j^JAf^nO 



(34) 



sin(7ri^)a„ 

where (i) is the Bessel function of the second kind in the standard notation and (t) 
is the Anger function |23j . 

Thus we expressed kernels K^(ri,r2) by combinations of S{ri — r2)l, (5(ri — r2)/3 
and G>o(t'i, T'2) operators, see Eq. (0). Operator 6o(t"i, ^2) of Eq. ([^ can be rewritten 
in a block form 

V T2(r-i,r2) -Ti{ri,r2) J 



where 



n— ^ ^ 



(36) 



71, n — 
•ii — y> 

Using Eqs. ([TT]) and ([20)1 we also can rewrite kernels A'^(ri,r2) in a block form 
where we denote (5 = 5{ri — r2)l. Then the kernels are 



A:*(ri,r2) 

/ (5(ri-r 

E 



2V^ 



n=0 



/ (5(ri-r2) (Da - iDi)A„_i^„_i - V^i:)2A„_i,„ 

(5(ri-r2) -V^L'2A„,„_i -(Dg + ii:)i)A„,„ 

(Dg + ii:>i)A„_i,„_i -V7^i:)2A„_i,„ (5(ri-r2) 

V -y^D2An,n-l {-D3+iDi)A„,n <5(ri - r2) 

Equation (jSH) together with Eqs. ([Ml) - (IMl), dSQ]) - (|32]) are the final expressions 
for the matrix elements of kernel operator. The kernels A'*(ri,r2) are given 
by combinations of two parts: translational-invariant gauge-independent parts and 
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the gauge-dependent phase factor e^^. Since x is not translational invariant, the 
kernels K^{ri,r2) are functions of independent variables ri and r2. 

The non-diagonal elements of matrix p9p are non-local in ri and r2 coordinates. 
In consequence, they change shapes of functions transformed with the use 
of K^{ri,r2), see Eq. pO)) . The transformed functions \^'^(ri) are more delocalized 
than the initial functions ^'(ri), since the non-diagonal elements of matrix (p9|) are 
characterized by the finite extent on the order of Ac or of magnetic length L. We 
describe this in the next section. 

The kernels K^{ri,r2) depend on two independent variables, i.e. on six 
coordinates. This case is more complicated than the field-free situation, where 
the transformation kernels depend on one independent variable R = ri — r2, see 
Refs. [niHn]- to slmpUfy our analysis, we set r2 = and denote ri = r. Using this 
simplification the matrix ([M)) can be written as 

( 5{r) f3(r)~zfi(r) ~f2(r) \ 

5{r) ^T2{r) -T3{r)-iT^{r 

f3(r)-Hzfi(r) -f2(r) 5{r) 

V -T2{r) -T^{r)+iTi{r) 5{r). J 

In the above equation instead of Di and A,„_„ integrals we defined another three 
functions 

ri(r) = \^^2 (41) 

n=0 

r2(r) ^ ^^^-"J ^2 E^o(«"ICI)LLi(p'), (42) 

^ n=l 



(40) 



r3(r) 



n— 



-+A(anC) 

z 



n— 

We introduced the notation — {x^ -f- y^) j (^l?) . In the second line of 
Eq. (j43l) we expressed the infinite sum over Laguerre polynomials by ^(p) function, 
see [Appendix B[ Functions ri(T) and r3(r) are defined similarly to ri(r), and r3(r) 
replacing an — )■ fln+i. The function r2(r) can be obtained from r2(T") replacing ri^o — > 
ro,i, see Eqs. ([27l) - ([28| . 

The summations over n in Eqs. (j4ip - (j43l) extend to infinity but in practice 
they are limited to finite numbers of terms by an exponential decrease of i4ro(a„|C|) 
and AiflnC) with n. However, for low magnetic fields there are large numbers of 
Laguerre polynomials involved in the summations so we truncate them for iiro(a„|^|) < 
10~^ and |y^(a„C)l < 10^^ or, alternatively, for n > 10^. This approximation allows 
us to perform calculations of Ti{r) for B > 10^ T, see below. For strong magnetic 
fields there is 77 — > 00, so the arguments a„C = {z / Xc) y^^+rfn of the Bessel or Anger 
functions in Eqs. (HT|) - P5|) are large for n ^ 0, and the terms including oqC do not 
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Figure 1. Calculated contour plot of ri(r), as given in Eq. II41I I. for the 
plane y = 0. Contour lines start for ri(7') = 0.05A^^ (outer lines) and end 
for ri(r) = lOAJ^. Coordinates x and z are in Ac units. For low magnetic fields 
the kernel is almost symmetric in x, y, z coordinates, while for large B the kernel 
in X and y directions is much more confined than in z direction. 



depend on the field. Then the summations over n reduce to single terms with n — 0. 
In this limit there is r2(r) ~ 0, while 



e 



-p^/2+ix 



ri(r)^ ttAL^ ^o(kl/Ac), (44) 



r3(r-) - + -^(^/Ac). (45) 



iS{p^) ie-P'/^+'>^ 
ttXcL'^z ttXcL"^ 

Thus both ri(r) and r3(r) have the extent L in the x — y directions and Ac in the z 
direction, so they are strongly anisotropic. For ri(r), r2('') and r3(T) functions there 
are no field- independent terms in the series (because oq — oi), so all these terms 
vanish in high magnetic fields and no widening occurs. 

In Figure[T]we show a contour plot of Fi (r) function for the y = plane calculated 
for various magnetic fields. For each of the four sections the contour lines start 
from Fi(r) = O.OSAj^ and they increase logarithmically to lOA^, ^. For all magnetic 
fields and all x and z the values of ri(r) are positive and they monotonically decrease 
to zero with growing x and z. For z = the expressions for Fi(r) diverge [see Eq. (PT|)]. 
but for X — and nonzero z they are large but finite. 

In Figure[T^ we show Fi (r) for a magnetic field B — 4.4x 10^ T, which corresponds 
to 1% of the Schwinger field of Bq = 4.4 x 10^ T. For such a field Fi(r) is, to a 
good approximation, spherically isotropic. The characteristic width of Fi(r) is on 
the order of Ac in all directions, as expected for the zero-field limit of the Dirac 
equation ^9j. Similar results are obtained for lower magnetic fields. At the Schwinger 
field (Figure [TJj) the term Fi(r) is still almost spherically isotropic but the widening 
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of ri(r) is larger than in the low-field limit. For still higher fields the terms ri(r) 
become anisotropic: their extent in the z direction changes slowly but they significantly 
shrink in the x — y directions. The behavior of ^2{f) and r3(r) functions is similar to 
that described for Fi (r), however F2(r) vanishes in the limits of i3 — >■ oo or i? — 0. For 
magnetic fields lower than the Schwinger field the behavior of Fj(r) terms is similar 
to the analogous Fj(r) terms. For higher fields all these terms go to zero. 



3. Second moments of transformation kernels 



Now we estimate quantitatively the spatial extent of the elements of kernel 
matrices K^{r,0) given in Eq. pO)) by calculating second normalized moments Wj 
of Fj(r) and Tj{r) functions (j=l,2,3). Here we follow Rose who used second 
moments to characterize non-locality of the transformation kernel for the Foldy- 
Wouthyusen transformation in absence of external fields. The second normalized 
moments are defined as 

/ r'^rj{r)d^r 
/r,(r)d3r 

The second normalized moments Wj of ^j{r) functions are defined similarly. The 
above integrals over F2(r), r3(r), F2(r) and F3(r) functions vanish due to symmetry. 
Below we analyze in more detail the second moments Wi and Wi of the ri(r) 
and ri(r) functions. Introducing auxiliary quantities 

Mt^^ f z''Ko{\z\aJ\c)dz, (47) 



where p = p/{-j2L), we obtain 



(48) 



2 kAP 

2 /„3 



Direct calculations give A^g — '^/o-n and A^f = A^/a„, where a„ are defined in 
Eq. ([33| . On setting in Eq. pS]): v = Q and t = p^, making the substitution xy/ (2L^) = 
sin(20)/(4L^) and integrating over d(j), we obtain 

/•oo 

Xg= / e-*/Vo(t/2)L°(t)dt - \/2P„(0), (50) 
where P„(a;) are the Legendre polynomials. Using the above results, we have 

f:A^gA^s=v2f:^, (51) 



n— 



Y^MlM^,^V2XlY.-^. (52) 
To calculate proceed in a similar way and obtain 



= 2L2 / e-*/Vo(i/2)tL°(i)dt. (53) 
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In order to use the result ([50)) we must express the product in terms of L°„(i) 

polynomials. Using the identities: (n + l)L^_^_^{t) - {2n + 1 - t)L° (t) + nL^^_^{t) = 
and L^{t)t = L^{t) - L^(t) we have 



or2 /-oo nr2 poc 

^ = / e^''^Jo{t/2)hl{t)tdt + / e-'/^Mt/2)hl{t)tdt 

f— ; an Jo o-o Jo 



n=0 n—1 



2i^V2<' Vp„(0) {---^) + VnP„(0) _ ^ _ ^ ) ^ , (54) 



I \0.n Q-n+l / — „ \CI71 0,n+l Q-n-l 

where we put l/a_i = and used the identities: P2n+i(0) — and nP„(0) = for n = 
0, 1. The sums over n can be calculated numerically either noting that for m — 2n there 

is Pm(0) — (—1/4)™ ( 1 or, alternatively, by a direct summation with the use of 



the generating functions for the Legendre polynomials, see [Appendix A Calculating 



the second moment of f i(r) one should substitute a„ a„+i in equations ([5T[) . ([52[) 
and ([54[). 

In the limit of high fields there is r/ ^ 1 and l/a„ ~ Sn^, so the second 
moment of ^i{r) in the x — y directions is Wi ~ 2L^ cx B~^, while the second 
moment of ri(r) in the z direction is Wf ~ A^, see Eq. ([44[) . For Wi there are 
no field-independent terms in the summations in Eqs. ([FT[) . (15^ and ([M|. and in 
high fields there is VV^ ~ Wf ~ 0. Thus no widening occurs. For low magnetic 
fields, when 77 ^ 1, we may approximate a~P ~ 1 ~ prfn. Since X]nPn(0) — 
(see [Appendix A I we have from Eq. ([5T[) that lims^o X^^^o -^o-^o = 1- This gives. 



upon using Eqs. ([49|), ([52|) and ([54[) . 

00 

lim ~ \/2A^ ^ P„(0)(1 - 3772^/2) ~ A^, (55) 

ra=0 

lim Wf ~ 2L2v^ <^ V P„(0)^ \ = 2\l, (56) 



The same results are obtained for Wf and moments. Thus the second moments 
of ri(r) and V\(r) functions in the low-field limit are 

lim Wi = lim Wi = 3A2. (57) 

This result reproduces our value obtained previously for MOT at _B = 0, see Ref. |10) . 

In Figure [2) we show Wi (thick line), dashed line) and Wf (dash-dotted line), 
as defined in Eqs. (HHl) . ([Fl]) . ([5^ and ([M[) . versus magnetic field B. For low magnetic 
fields WC tends to 2A^ and remains nearly constant up to the Schwinger field, where 
a flat maximum occurs. For larger magnetic fields starts to decrease and finally 
decreases as . The second moment remains nearly constant, Wf ~ A^, within 
the whole range of magnetic fields except in the vicinity of the Schwinger field. The 
total second moment Wi = Wf -t-Wf of the ri(r) function, exhibits for an increasing 
magnetic field a transition from the 3D value Wi = 3A^ to the ID value Wi — A^. 
For low magnetic fields and Wf moments are close to and . However, for 
magnetic fields higher than the Schwinger field both of them tend to zero exponentially 
with B. It should be reminded that the kernels given by Eq. P0[) are composed of 
various elements. All of them are nonlocal (with the exception of the delta functions 
on the diagonal), but only ri(r) and ri(r) have non- vanishing second moments. It is 
the second moment of ri(r) that is shown in Figure [2) 
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Figure 2. Calculated normalized second moment of ri(7'), as given in 
Eq. I|49| l. versus magnetic field B. Solid line: total second moment Wi, dashed 
line: Wj , dot-dashed line: Wj . Second moment Wi exhibits a transition from 3D 
case Wi ~ 3A|? at low fields to ID case Wi ~ at high fields. 



For a non-relativistic electron in a magnetic field there is no widening of the 
transformed wave function. We retrieve this limit by setting mc^ — > oo, which 
gives: — > 0, — > and a,i = 1 for all n. Then the sums in Eq. vanish, 
the summation in Eq. (jSip gives unity, and the summation in Eq. (j52p gives — > 0. 
Therefore = 0, VVf = 0, VVf = and no widening occurs. 



4. Transformed functions 



In this section we consider properties of transformed functions using as an example 
a Gaussian function of finite width. We take the initial function 'i'{r) in the 
form: ^'(r) = /(r)(0, 1, 0, 0)^, in which 

/M = -Vid^' (58) 

and d characterizes function's width. The transformed function is: ^'{r) = V'i>{r) = 
/(r)/v/2(0, 1,0, i)^- Then we have l^-'*) = which gives [see Eq. ^] 



V2 



/ 



± 



2 ^ 

n=0 




firi 


-iy^D2{zi - Z2)A„_l,„(pi, P2) 
{£'1(21 - Z2) - iDsizi - Z2)}An,„{pi,P2 

-y/nD2{zi - Z2)A„_i,„(pi, P2) 

\ {iDi{zi- Z2) - D3{zi- Z2)}An,„{pi,P2) J 



/(r2)d'r2,(59) 
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where A,„ „ and Di are defined in Eqs. - and ([50]) - ([5^ . As seen from 
Eq. ([5^ . positive or negative energy functions ^'^(ri) consist of two parts: the 
function ^''(ri) and the nonlocal parts arising from the integration of ^''(ri) over 
appropriate elements of r2). The summation over n and the threefold integrals 

over r2 can be performed numerically. 

Now we turn to the high-field approximation, in which L ^ A, 
the summation over n in Eq. 
transformed function is 



59|) reduces to the single term with n 



In this case 
= and the 



*'±(r) 



/2C 
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/ 






f{r) 


± - 






Fz{z)Fp{p) 
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iF,{z)Fp{p) j 



where C is a normalization constant and 

/2d /"^ (1 - acfc^)e-'=''*'/2+''=^ 



27r3/4 



The presence of Gaussian term ensures the convergence of integration. Further 

1 



(60) 



(61) 



Fp{p) 



2V27rd 



P2 



V(F+2(Pl2 + 2l4 



-a^x —a-yV +«o xy 



where 



2(d4 + 2d2L2 ^ 



2L4)' 



6^ = 



2L2(d4_ 



2d^L^ + 2L4) ^ 



d4 + 2^2^2 + 2X4' 



For narrow initial functions: d L, there is 



(62) 

(63) 
(64) 
(65) 

1/(4L2) and b"^ ~ l/(2£2), 



so that Fp{p) has a spatial extent on the order of L/2, similar to the case analyzed 
before, see Eqs. (|311) and P5|) . In the opposite limit of wide initial functions: d ^ L, 
we have ~ &2 ^ 1/(2^^) and Oy ~ 1/{2L'^). The differences between Ox and Cj, 
result from the asymmetry of the Landau gauge. To obtain the same results in another 
gauge, one would have to modify the initial functions accordingly, see Refs. [26l [27] . 
To demonstrate a widening of the functions 'i''^{r) we calculate their variances 



= j 'i!'^{r)^rH'^{r)d^r j J \^'^{r)\^d^r , 



(66) 



which are different from the second moments calculated above, see Eq. (|46p. The 
integrals in Eq. can be calculated analytically but we do not quote the long 
resulting expressions. 

The plots of for various magnetic fields and initial function widths are given 
in Figure [3l The solid lines represent calculated variances of the transformed function. 
The variances Zc = Sd'^ /2 of the initial functions are indicated on the right ordinate. 
For the initial function of the width d the variances of the transformed functions are 
larger than Zq, so the transformed functions are wider than the initial one. For the 
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Figure 3. Calculated variances of transformed Gaussian function ^'^(r), as 
given in Eq. Il60|l . versus magnetic field B/Bq in the high-field approximation. 
Solid lines: variances corresponding to the initial functions having different 
widths d. Variances of the initial Gaussian functions for d = Xc and d = 0.5Ac 
are indicated. 

initial functions having d ~ Xc the variances reach their limit Zq for magnetic 
fields on the order of = 10^^ T. For smaller d, the variances of transformed functions 
are much larger than Zq and they reach the limit Zq for very large fields. All results 
presented in Figure [3] are valid for magnetic fields larger than the Schwinger field, for 
which L ^ Ac- The asymmetry between positive and negative variances Z^ shown 
in Figure [3] arises from the asymmetry in the initial wave function having only one 
nonzero component. Selecting other nonzero component or more than one nonzero 
components, one would obtain different numerical values of Z^, but their general 
properties would be similar to those presented in Figure [3l 

5. Tv^ro-dimensional case 

Recently, Lamata et al. j21| proposed new topological effects for the two-dimensional 
Dirac equation in a magnetic field. They also indicated a possibility of simulation 
and reconstruction of the Wigner function for the lowest Landau levels with the use 
of trapped ions. In this connection it is of interest to calculate the non-locality of the 
transformation kernel in two dimensions (2D). The matrix elements of 2D kernels can 
be obtained from the three dimensional results given above by setting e'*^^^ — )■ 6{kz) 
in all integrals including kz and by a proper adjustment of normalization constants. 
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Figure 4. Gaugc-indcpcndent part of kernel's element Gi(p) in two dimensions, 
as given in Eq. JTO}, calculated for four values of magnetic field. For low fields the 
kernel non-locality is on the order of Ac and it does not depend on the field. For 
higher fields the non-locality slightly increases and for very large fields it decreases 
again. 



The 2D counterpart of Eq. gOl) is 

^/3± i 



/ S{p) 

6{p) 

+iGi{p) -G2{p) 

\ -G2{p) +iGiip) 



-iG^{p) 

~G2(p) 

5{P) 




-iG,{p) 


<5(P). / 



(67) 



In the above expression we defined 
Gi(p) 



G2{P) = 



L2 



E 



(69) 



using the same notation as in Eq. (|40p . Functions Gi{p) and G2{p) are defined 



similarly to Gi{p) and G2{p) replacing 



rl^{n + 1), see Eqs. ^ 



Note that in the 2D case there is no counterpart of the r3(r) term. 
Functions Gi{p), G2(p), Gi{p) and G2{p) are products of gauge-independent parts 
and the gauge- dependent factor e'^. The summations over n in Eqs. (|68t - (|69| 
are slowly convergent. However, these sums can be transformed to fast converging 
integrals that can be calculated numerically. For Gi{p) one has 
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where t = exp(— 77^(7^). Details of the derivation are given in [Appendix A[ The 
remaining sums can be transformed into integrals in a similar way. 

In Figure |3] we plot the gauge- independent part of Gi{p) for four values of 
magnetic field. For low fields the spatial extent of the kernel is on the order of Ac 
and it does not depend on the field. For the Schwinger field (dash-dotted line) the 
non- locality of Gi{p) differs only slightly from its low-field value (solid line). For 
larger fields the non-locality of Gi{p) increases (dashed line) and for very large fields 
it decreases again (dotted line). The non- locality of Gi{p) is on the order of the 
magnetic length L which in this regime is much smaller than Ac. For large fields in 
three dimensions the non-locality of the kernels tends to zero in the x and y directions, 
but it remains finite in the z direction. For large fields in 2D the non-locality of the 
kernels goes to zero. Therefore, for large fields in 2D there is no widening of the 
transformed functions. This conclusion may be important for simulations of the Dirac 
equation by trapped ions 12 Ij . where arbitrary large values of a simulated magnetic 
field may be achieved. 

6. Discussion 

As follows from our considerations and, in particular, from the above figures, the 
magnetic field effects for relativistic electrons in a vacuum become appreciable at 
gigantic field intensities, comparable to the Schwinger field Bq — 4.4 x 10^ T. 
However, one can overcome this problem in two ways. The first is provided by 
narrow-gap semiconductors, whose energy band structure can be described by Dirac- 
type Hamiltonians with considerably more favorable parameters, see Ref. [151 • For 
example, in InSb the energy gap, corresponding to 2mc^ in a vacuum, is Eg ~ 0.23 
eV and the effective mass m*, corresponding to rest electron mass m in a vacuum, 
is m* ~ 0.014rn. In consequence, an effective Schwinger field in InSb, given by the 
condition heBQ/m* — Eg/2, is Bq ~ 14 T. Such magnetic fields are easily available 
in laboratories, so the field effects described above can be readily investigated in 
narrow-gap semiconductors. One should, however, bear in mind that in experiments 
with crystalline solids one deals with additional effects, most notably with electron 
scattering leading to broadening of Landau levels, which were not taken into account 
in our considerations. 

The second way to lower effective magnetic fields is to use simulations of the 
Dirac equation employing trapped ions or cold atoms interacting with laser radiation, 
see e.g. Refs. |T71 [THl [H] and review [7]. Such simulations may also include the 
presence of an external magnetic fields [2j [21] . In this type of experiments an effective 
simulated field can be tailored to suit observation requirements. For example, in a 
simulation of the Dirac equation proposed in Ref. [2] , the effective simulated magnetic 
field corresponds to the ratio (fieB / ra) l(2mc^) = 16.65 for realistic trap parameters. 
However, in simulation- type experiments there is no need to apply real magnetic fields 
apart from the fields necessary to create ion traps. 

As to our results, the general rule is that, similarly to the non-locality of 
energy separating transformations for i? = 0, the characteristic quantity is the 
Compton wavelength Ac = h/mc, see Ref. [lO]. In the presence of a magnetic 
field, typical anisotropy appears between the z direction along the magnetic field 
and the transverse x — y plane. This is seen in the broadening of the transformation 
kernel illustrated in Figure [1] The anisotropy appears also in the second moments of 
functional transformation kernels shown in Figure [21 This result can be considered to 
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be a direct extension of that of Rose [5], who characterized quantitatively the non- 
locahty of the Foldy-Wouthuysen transformation for B = 0. The low field limit Wi = 
3A^ agrees with our previous calculation for MOT at B = 0. The corresponding result 
for FWT is Wi = (3/4) and one can say that the FW transformation is somewhat 
more "compact" than the MO transformation. Interestingly, at high fields Wi is 
suppressed, which results in Wi — . This means that at high B the quantization of 
the electron motion in the x—y plane gives in consequence an effective one-dimensional 
motion along the magnetic field. The parallel second moment Wf depends somewhat 
on B because in the relativistic mechanics there is no exact separation between parallel 
and transverse directions of magnetic motion. Finally, Figure [3] indicates two effects. 
First, the broadening of transformed functions is stronger for narrow initial widths. 
Second, with an increasing magnetic field the broadening of transformed functions 
diminishes, reaching asymptotically the non-broadened values at very high fields. One 
should bear in mind that the results shown in Figure [3] were obtained in the high field 
approximation, i.e. B / Bq > 1. 

The results obtained in this paper refer to a non-locality of the energy-separating 
transformation for the Dirac Hamiltonian in the presence of a uniform magnetic 
field. However, under two assumptions given below it is possible to generalize our 
approach to the case of an inhomogeneous magnetic field. First, there should exist 
an exact transformation separating the positive and negative energy states. Second, 
it is practical to have analytical solutions of the Dirac equation for electrons in an 
inhomogeneous magnetic field that allow one to find an expansion of the kernel into 
eigenvectors in a way analogous to that given in Eq. ([T2l) . More specifically, the 
existence of transformation separating positive and negative energy states can be 
simplified to an assumption that the operator Ti ^ [see Eq. Q], describing the square 
of the Dirac Hamiltonian in an inhomogeneous magnetic field, is diagonal. This 
criterion can be satisfied for example by a vector potential having the form A = 
[Ax{x,y), Ay{x,y),0] with arbitrary functions Ax{x,y) and Ay{x,y). Then the 
operator separating the positive and negative energies is analogous to U given in 
Eq. ([5]) but with a suitably changed vector potential. Physically, the existence of such 
an operator is limited to situations in which there is no creation of electron- hole pairs. 
Therefore, such an operator does not exist for Dirac particle in an electric field or in 
some configurations of the magnetic field that mix positive and negative energy states. 
As to the second requirement for the applicability of the present approach, there are 
several inhomogeneous magnetic fields known in the literature for which analytical 
solutions exist. As an example, in three dimensions the Dirac equation has solutions 
for vector potentials A = [(B/a) tanh(ay), 0, 0] or A = [{B/a)[l — exp{ay)],0,0], 
where is a characteristic length of the potential |28j . In two dimensions there 
exist also solutions for magnetic field B = /3 sec^ {ay) . The above vector potentials 
have the form A — [Ax{x, y),Ay{x, y),0], so that the operator H ^ is diagonal and one 
can calculate the non-locality of the transformation kernel in a way analogous to that 
presented above. We expect the results to be similar to those presented in our work 
with properly modified magnetic length. 

Our calculations of the non-locality generated by FWT and MOT at S = 0, see 
Ref. [lOj . indicate that different "separating" transformations have similar properties 
and small differences between them concern mostly the extent of non-localities. 
In consequence, we believe that the properties of MOT investigated in this work 
characterize other separating transformations for the Dirac electron in the presence of 
a magnetic field. 



Non-locality of transformations in a magnetic field 



18 



7. Summary 

We investigate non-locality of transformations separating positive and negative energy 
states of the Dirac electrons in the presence of a magnetic field, using as an example the 
Moss-Okninski transformation. Functional kernels generated by the transformation 
are described and their second moments in the coordinate space are calculated to 
determine quantitatively the nonlocal features. The latter are characterized by the 
Compton wavelength Ac = h/mc. The behavior of kernels is different along the 
magnetic field and in the plane transverse to it. Dependence of non- locality on the 
intensity of magnetic field is studied. It is demonstrated that the transformed functions 
are broadened. Finally, it is indicated that a practical way to study Dirac electrons in 
the presence of a magnetic field is to simulate their properties using either narrow-gap 
semiconductors or trapped ions interacting with the laser radiation. 



Appendix A. 

Calculating integrals Am,n in Eqs- ((26|) - p8| we use the identity 

/oo 
H„,((7 + a)H„((j + b)e-'''dq = 2"V^m!6"'L;^' (-2a&), (A.l) 
-oo 

where n' = n — m and m < n. 

A convenient way to calculate sums in Eq. ()54|) is to use the generating 
functions for the Legendre polynomials: X^i^o — ~ ^te -I- 

and X^i^o ~ ~ ~ 2te -1- t^)'^/^. In order to apply the generating 

functions we first make the transformation l/a„ = 2/y/TT exp[—q^ — q^ri^n]dq, next 
we set t = e~'^ , then change the order of summation and integration. As a result 
one obtains 

^ Pn(0) - 2 e-o'dq 

nP„(0) r e'^'-^^'^'dq 

For vanishing 77 (weak magnetic fields) the above integrals tend to l/\/2 and — 1/\/8, 
respectively, while for 77 1 (strong magnetic fields) they tend to ±1. 

The same method can be used to calculate a sum of the type Xn-Pn(0)/(1 + 
Ty^n)'^/^, employing the identity 1/(1 -I- r/^n)'^/^ _ —2[d/db)\b-\-Tfn]~'^/'^ taken at = 1. 
Proceeding the same way as for Eq. \A.2\ we get 

P„(0) ^ J_ r g^e-^^dg 



which applies to Eq. (|52l). 

Consider the sum in Eq. (|68| for the 2D case. Using the identity: 1 / -^/l + rfn = 
(1/0F) e-«'-™''«'dg we obtain 

s = f; = 4 i: /" ^sf)'-'' U-''")' 
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Changing the order of summation and integration, and using the identity 

flKi^r = ^^3^exp (A.6) 

with t — exp(— 77^(7^) one obtains Eq. ([70)1 . The presence of a Gaussian term in 
Eq. (|A.5I) ensures fast convergence of the integral for large q. For smaU q the 
factor exp[te/(t — 1)] tends to zero as exp(— l/q^) and the integrand in Eq. (155)) 
remains finite for all values of q. 



Appendix B. 

We prove a sum rule for the A„^„ integrals. Let pi and p2 be position vectors in 
the two-dimensional space and {|q)} the complete set of eigenstates of a hermitian 
operator. Then 

^(Pl - P2) = (P1|P2> = 5](pi|q>(q|p2> = 5] Vq(Pl)V'q(P2)^ (B.l) 

q q 

where V'q(p) = (pIi)- Taking the eigenstates |q) = \n,kx) of the Hamiltonian for a 
non-relativistic electron in a magnetic field we obtain 



/ on ^ rr H„(a)H„(e2)dfc.-<5(pi-p2), (B.2) 

where p = {x,y). Performing in Eq. (jB.2p the integration over kz and using the 
notation from Eq. (|25|) we have 



E 

ra=0 



n.n 



5{pi - P2). (B.3) 



Identity (jB.3[) was used in Eq. for the calculation of infinite sums over the Laguerre 
polynomials appearing in the l/z term. 
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